FIRST ATTEMPT AT MEASURING THE CMB CROSS-POLARIZATION 
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We compute upper limits on CMB cross-polarization by cross-correlating the PIQUE and Saskatoon 
experiments. We also discuss theoretical and practical issues relevant to measuring cross-polarization 
and illustrate them with simulations of the upcoming BOOMERanG 2002 experiment. We present 
a method that separates all six polarization power spectra (TT, EE, BB, TE, TB, EB) without any 
other "leakage" than the familiar EE-BB mixing caused by incomplete sky coverage. Since E and 
B get mixed, one might expect leakage between TE and TB, between EE and EB and between BB 
and EB — our method eliminates this by preserving the parity symmetry under which TB and EB 
are odd and the other four power spectra are even. 
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I. INTRODUCTION 

Although not yet detected, it is theoretically expected 
that the Cosmic Microwave Background (CMB) is polar- 
ized. CMB polarization should be induced via Thomson 
scattering which occurs either at decoupling or during 
rcionization. The level of this polarization is linked to the 
local quadrupole anisotropy of the incident radiation on 
the scattering electrons, and it is expected to be of order 
1%-10% of the amplitude of the temperature anisotropics 
depending on the angular scale (see jl]|2| and references 
therein). 

CMB polarization is important for two reasons: hrst, 
polarization measurements can substantially improve the 
accuracy with which cosmological parameters are mea- 
sured by breaking the degeneracy between certain pa- 
rameter combinations; second, it offers an independent 
test of the basic assumptions that underly the standard 
cosmological model. 

Since the polarized CMB signal is so small, it is quite 
likely that its first detection will be an indirect statisti- 
cal one, from its predicted correlation with an unpolar- 
ized CMB map. Since such cross-correlations involve one 
rather than two powers of the (small) polarization frac- 
tion, they will be measured with better signal-to-noise 
than polarization autocorrelations. Indeed, it has been 
shown HQ that for almost all cosmological parameters, 
the polarization capabilities of upcoming high-precision 
CMB experiments such as MAP and Planck add infor- 
mation mainly through the cross-polarization signal, the 
only exceptions being the reionization and gravity wave 
parameters. 

The goal of the present paper is to present a detailed 
study of the cross-polarization from a practical point of 
view, connectingreal- world data to physical models. In 
section Section nil we argue that the dimensionless cor- 



relation coefficient rg is a more meaningful quantity to 
discuss than the cross power spectrum Cf , and illus- 
trate how it depends of various cosmological parameters 
for standard adiabatic models. 



In Section III, we compute the strongest observational 
constraints to date on rg by cross-correlating the polar- 
ized Princeton IQU Experiment (PIQUE) S with the 
unpolarized Saskatoon (SK) data set. We do this 
using the formalism presented in H , which takes into ac- 
count real-world issues such as incomplete sky coverage 
and correlated noise. Finally, in Section IV we assess the 



prospects for measuring rg in the near future by analyz- 
ing a simulated version of the upcoming BOOMERanG 
2002 experiment. 



II. POLARIZATION PHENOMENOLOGY 

Whereas most astronomers use the Stokes parameters 
Q and U to describe polarization measurements, the 
CMB community uses two scalar fields E and B that 
are independent of how the coordinate system is ori- 
ented, and are related to the tensor field (Q, U) by a 
non-local transformation p| . Scalar CMB fluctua- 
tions have been shown to generate only ^-fluctuations, 
whereas gravity waves, CMB lensing and foregrounds 
generate both E and B. This formalism has been ap- 
plied to real data by the PIQUE § and POLAR Q 
teams. 
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A. The six power spectra 



Since CMB measurements can be decomposed^ into 
three maps (T,E,B), where T denotes the unpolarizcd 
component, there are a total of 6 angular power spectra 
that can be measured. We will denote these Cf, Cf, 
Cf, Cf, Cf and Cf , corresponding to TT, EE, BB, 
TE, TB and EB correlations^, respectively. By parity, 
CY = Cf = for scalar CMB fluctuations, but it is 
nonetheless worthwhile to measure these power spectra 
as probes of both exotic physics |l^-[l8|] and foreground 
contamination. Cf = for scalar CMB fluctuations to 
first order in perturbation theory ||^- [lfj| , |l9| — secondary 
effects such as gravitational lensing can create B polari- 
zation even if there are only density perturbations present 
p0| . The remaining three power spectra are plotted in 
Figure [I] (top) for the "concordance" model of (that 
of ^| is very similar), showing that Cf is typically a 
couple of orders of magnitude below Cj on small scales 
and approaches zero on the very largest scales (in the 
absence of reionization). 




FIG. 2. The region of the parameter space (r 



) allowed 



by the generalized Schwarz inequality. The cube extends from —1 
to 1 in all three dimensions. The four corners (1, 1, 1), (1, —1, —1), 
(—1, 1, —1) and (—1, —1, 1) correspond to T, E and B being per- 
fectly correlated and anticorrelated. 
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FIG. 1. The T, X and E power spectra are shown (top panel) 
for the concordance model of J2l| . The TE correlation coefficient 
rf (bottom panel) is the ratio between Cf and the geometric 
mean of Cf and Cf , so in the top panel, \rf | is just the distance 
between TE and the geometric average of the T and E curves 
(dashed curve). 



* Since the transformation between (Q, U) and (E, B) is non- 
local, the ^/^-decomposition is straightforward only for the 
case of complete sky coverage. Methods for proceeding in 
practice for the real-world case of incomplete sky coverage 
are discussed in j^,^o||l^-|l^] , and we will use such a method 
for our real-world calculations below. 

t From here on, we adopt TT = T, EE = E, BB = B, 
TE = X,TB = Y, EB = Z, 



B. Covariance versus correlation 



The cross-power spectrum Cf is not well suited for 
such a logarithmic plot, since it is negative for about half 
of all ^-values. A more convenient quantity is the di- 
mensionless correlation coefficient rf plotted in Figure ^ 
(lower panel), defined as 
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(i) 



since the Schwarz inequality restricts it to lie in the range 

(2) 



1< rf < 1. 



These dimensionless correlations/anticorrelations are 
seen to be quite strong in the sense of being near these 
limiting values on many scales. For instance, if one were 
to smooth the T- and i?-maps to contain only large an- 
gular scales I < 10 where rf ~ 0.9, they would be so 
strongly correlated that most hot and cold spots would 
tend to line upf[ Conversely, if one were to band-pass 
filter the two maps on scales 100 < I < 200 where 
rf < —0.6, hot spots in the T-map would tend to line 
up with cold spots in the iJ-map. 

More generally, let us also define rf = Cf /(CfCf) 1 / 2 , 
rf = Cf /{CfCf ) 1 ' 2 . Expanding the T, E jmd B 
maps in spherical harmonics with coefficients a 
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* Such a correlation, however, would be very difficult to de- 
tect given how small E is expected to be in that l?-range 
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and af m , the three cross power spectra are by defini- 
tion the covariances Cf = (aj^afj, Cj = (aj^afj, 
Cf = (afm a f m )> so rX ' rY an d rZ are simply the corre- 
lation coefficients between aj m , af m and af m . 

What values are these three correlation coefficients al- 
lowed to take? They are real-valued just as T, E and 
B, which is most easily seen using real- valued spherical 
harmonics. The dimensionless correlation matrix corre- 
sponding to the vector (aj m , af m , af m ) is 

fl rfrj\ 
R= rf 1 rf , (3) 
VI rf 1 J 

and by virtue of being a correlation matrix, it cannot 
have any negative eigenvalues. This implies not only that 
\ r f\ < 1, \ r J\ < 1 an d \rf\ < 1, but also that the deter- 
minant of R must be non-negative, i.e., that 

{rff + {rjf + {rff -2rfrjr z l <1 (4) 

for all I. This allowed region in the 3-dimensional space 
(r x ,r Y ,r z ) is plotted in Figure |[ and is seen to re- 
semble a deformed tetrahedron. If two coefficients are 
strongly correlated with each other, then they must both 
have roughly the same correlations with the third, ap- 
proaching the limiting case (r x ,r Y ,r z ) = (1,1,1) (up- 
per right corner). The other three corners correspond 
to the allowed possibilities (1,-1,-1), (—1,1,-1) and 
(-1,-1,1)- 

From here on we use as shorthand for r x . 

C. Cosmological parameter dependence of 
polarization power spectra 

A detailed review of how CMB polarization reflects 
underlying physical processes in given in . In this sub- 
section, we briefly review this topic from a more phe- 
nomenological point of view (see also [@), focusing on 
how different cosmological parameters affect various fea- 
tures in the E, B and X power spectra and aiming to 
familiarize the reader with, in particular, the correlation 
spectrum m and its cosmology dependence. For more de- 
tails, the reader is refererred to the polarization movies 
at www.hep.upenn.edu/^ angelica/polarization. html. 

Let us consider adiabatic inflationary models specified 
by the following 10 parameters: the reionization opti- 
cal depth r, the primordial amplitudes A s , A t and tilts 
n 8 , n t of scalar and tensor fluctuations, and five param- 
eters specifying the cosmic matter budget. The various 
contributions f2j to critical density are for curvature f2k, 
vacuum energy Q\, cold dark matter O c d m , hot dark 
matter (neutrinos) f2„ and baryons fib- The quantities 
= h 2 £lb and cjdm = ft 2 fidm correspond to the physical 
densities of baryons and total (cold + hot) dark matter 
(f2dm = fi c dm + ^i/), and f v = fi^/fidm is the fraction of 
the dark matter that is hot. The baseline values of the 



parameters here and in the movies are for the concor- 
dance model of [SJ, r = £lk = A t — f u = 0, Qa = 0.66, 
^dm = 0.12, Wb = 0.02, n s — 0.93, providing a good fit 
to current data from the CMB, galaxy and Lyman Alpha 
Forest clustering and Big Bang nucleosynthesis. 

1. Polarized versus unpolarized 

If recombination were instantaneous, there would be 
no polarization at all! 

Both the E and the T power spectra carry information 
about the z > 10 3 pre-recombination epoch in the form 
of acoustic oscillations. From a practical point of view, 
there are two obvious differences between the E and T 
power spectra as illustrated by Figure |l|: 

• The E power is smaller since the polarization per- 
centage is small, making measurements more chal- 
lenging. This is because polarization is only gener- 
ated when locally anisotropic radiation scatters off 
of free electrons, and this only occurs during the 
brief period when recombination is taking place: 
before recombination, radiation is quite isotropic 
and after recombination there is almost no scatter- 
ing. 

• Aside from reionization effects, the E power ap- 
proaches zero on scales much larger than those of 
the first acoustic peak. This is because the pola- 
rization anisotropies are only generated on scales 
of order the mean free path at recombination and 
below. 

As detailed below, changing the cosmological parame- 
ters affects the polarized and unpolarized power spectra 
rather similarly except for the cases of reionization and 
gravity waves. All power spectra were computed with 
the CMBfast software [§5|. 

2. Reionization 

Reionization at redshift introduces a new scale 
~ 20(z*/10) 1 / 2 corresponding to the horizon size at 
the time. Primary (from z > 1000) fluctuations &Tg on 
scales £ 3> get suppressed by a factor e r and new series 
of peaks[|] are generated starting at the scale . Figure || 
illustrates that although these new peaks are almost un- 
detectable in T, drowning in sample variance from the 
unpolarized Sachs- Wolfe effect, the are clearly visible in 



s These new peaks are caused not by acoustic oscillations, but 
by a projection effect: they are peaks in the Bessel function 
that accounts for free streaming, converting local monopoles 
at recombination to local quadrupoles at reionization. 
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FIG. 3. How the reionization optical depth r affects the T and 
E power spectra (top) and the TE correlation rg (bottom). Solid, 
dashed and dotted curves correspond to for t=0, 0.2 and 0.4, re- 
spectively. 
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FIG. 4. How the gravity wave amplitude At affects the T, E 
and B power spectra (top) and the TE correlation ri (bottom). 
Solid, dashed and dotted curves correspond to for At=0, 0.2 and 
0.4, respectively. 



E since the Sachs- Wolfe nuisance is unpolarized and ab- 
sent. The models in Figure [| have abrupt reionization 
giving t oc z*^ 2 , so higher z„ is seen to shift the new 
peaks both up and to the right. 

On small scales, reionization leaves the correlation rg 
unchanged since Cj and Cf are merely rescaled. On 
very large scales, rg drops since the new polarized sig- 
nal is uncorrelated with the old unpolarized Sachs- Wolfe 
signal. On intermediate scales £ > I* ~ 20, oscillatory 
correlation behavior is revealed for the new peaks. 

For more details about CMB polarization and reion- 
ization see 26 1 . 



comes mainly from E, not T: since the tensor polariza- 
tion has a redder slope than the scalar polarization, it can 
dominate E at low £ even while remaining subdominant 
in T. 

The amplitudes A s , A t and tilts n s , n t of primordial 
scalar and tensor fluctuations simply change the ampli- 
tudes and slopes of the various power spectra: B is con- 
trolled by (A t , n t ) alone, whereas T and E are affected by 
(A s ,n s ) and (A t ,nt) in combination. Note that if there 
are no gravity waves (A t — 0), then these amplitudes 
and tilts cancel out, leaving the correlation spectrum rg 
independent of both A a and (apart from aliasing effects) 

71.,. 



3. Primordial perturbations 

As seen in Figure |], gravity waves (a.k.a. tensor fluc- 
tuations) contribute only to fairly large angular scales, 
producing E and B polarization. Just as for the reion- 
ization case, unpolarized fluctuations are also produced 
but are difficult to detect since they get swamped by the 
Sachs- Wolfe effect. As has been frequently pointed out 
in the literature, no other physical effects (except CMB 
lensing and foregrounds) should produce B polarization, 
potentially making this a smoking gun signal of gravity 
waves. 

Adding a small gravity wave component is seen to sup- 
press the correlation rg in Figure ^, since this component 
is uncorrelated with the dominant signal that was there 
previously. Indeed, this large-scale correlation suppres- 
sion may prove to be a smoking gun signature of gravity 
waves that is easier to observe in practice than the oft- 
discussed £?-signal. This Ti?-correlation suppression 



4- Spacetime geometry 

Just as A s and n s , the spacetime geometry parameters 
Ok and affect the polarized and unpolarized power 
spectra in similar ways. f2k and J7a were completely ir- 
relevant at z > 10 3 , when the acoustic oscillations were 
created, since Ok ~ Oa ~ at that time regardless of 
their present values. As is well-known and illustrated by 
the above-mentioned movies, the acoustic peak features 
are therefore independent of these parameters, merely 
shifting sideways on a logarithmic plot as geometric ef- 
fects magnify/shrink the scale of primordial fluctuation 
patterns. On large scales, the late ISW effect creates ad- 
ditional power that is completely unpolarized, since this 
is a pure gravity effect involving no Thomson scattering. 
Since the ISW effect is uncorrelated with the primary 
large-scale fluctuations, it suppresses rg on large scales 
as Ok or Oa shift away from zero as seen in Figure 
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FIG. 5. How the spatial curvature affects the T and E power 
spectra (top) and the TE correlation rg (bottom). Solid, dashed 
and dotted curves correspond to = —0.1 (closed model), (flat 
model) and 0.1 (open model), respectively. Apart from the familiar 
horizontal shift, the late ISW effect is seen to suppress the correla- 
tion ri on large scales. 
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FIG. 6. How the baryon fraction Qj,/(Qf, + Q^) affects the T 
and E power spectra (top) and the TE correlation ri (bottom). 
Solid, dashed and dotted curves correspond to baryon fractions of 
0.01, 0.02 and 0.08, respectively. 



A. Data 



5. Matter budget 

The Primordial CMB signal is dominated by fluctu- 
ations in the gravitational potential and the density at 
z ~ 10 3 , whereas the i?-signal is dominated by pecu- 
liar velocities on the last scattering surface. This is why 
the .E-spectrum is seen to be out of phase with the T- 
spectrum, peaks in one matching troughs in the other. 
This also explains why increasing the baryon fraction 
fb = f2b/(^b + ^dm) as in Figure || lowers the polarized 
peaks, in contrast to the boosting of odd peaks for the 
unpolarized case: more baryons lower the sound speed in 
the photon-baryon plasma, producing lower velocities. 

The remaining matter budget parameters, the cold and 
hot dark matter densities, effect E polarization in much 
the same way T. Increasing the dark matter density 
/i 2 fidm shifts the peaks down and to the left. The T, E 
and r£ power spectra change hardly at all when changing 
the fraction Qu/Qdm of the dark matter that is hot. This 
is because the neutrinos were already quite cold (nonrel- 
ativistic) at the time the CMB fluctuations are formed. 



III. CASE STUDY I: PIQUE 



PIQUE was a CMB polarization experiment on the 
roof of the physics building at Princeton University. 
It used a single 90 GHz correlation polarimeter with 
FWHM angular resolution of 0°.235, and observed a sin- 
gle Stokes parameter Q in a ring of radius of 1° around 
the North Celestial Pole (NCP) [§ (see Figure §). 

During one day, the telescope was able to observe on 
the order of 20 independent points on this ring, chopping 
slowly (every few seconds) between two points separated 
by 90° along this circle. The polarized sky signals de- 
tected at these two points had opposite signs (±Q) and 
were six hours out of phasep*]. This means that the ex- 
periment did not measure individual Q-values, but sums 
of two. To simplify subsequent calculations, we elimi- 
nated this complication using the deconvolution method 
described in appendix D of p4| to recover a filtered ver- 
sion of the Q"inap. Specifically, the ring was pixelized 
into 144 angular bins, and we denoted the corresponding 
Q- values Q\, Q144. Let us focus on four pixels form- 
ing a perfect square in the sky, say pixels 1, 37, 73 and 
109, and group them into a vector x. Because of the 
chopping, PIQUE measured not x but the linear trans- 
formation Ax, where 



We now turn to the issue of measuring cross- 
polarization in practice, including issues of methodology, 
window functions and leakage. We analyze an existing 
data set in this section, then turn to simulations of up- 
coming data in Section IV . 



See ]5| and dicke.princeton.edu for more details. 
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FIG. 7. The deconvolved PIQUE data consists of the Stokes 
Q-parameter as a function of position along a 1° radius circle 
around the North Celestial Pole. One quarter of the modes are 
projected out as described in the text. 
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The matrix A is singular, with a vanishing eigenvalue 
corresponding to the vector (1,-1,1,-1). Our decon- 
volution therefore sets this particular mode to zero, and 
encodes our lack of information about its value as near- 
infinite noise for this mode in the map noise covariancc 
matrix. Since there are a total of 36 such pixel quadru- 
plets, our deconvolved data set, which is plotted in Fig- 
ure [7J, thus has 36 such unmeasured modes. 

Because PIQUE was insensitive to the unpolarized 
CMB component, we cross-correlated the Q-data from 
PIQUE with T-data from the SK map Q, deconvolved 
and pixelized as described in |24|. To take advantage of 
cross-polarization information from spatially separated 
pixels, we used SK pixels not merely from the PIQUE 
circle, but from a filled disk of radius 3° around the NCP, 
a total of 288 0.31° x 0.31° pixels. We found that further 
increasing the size of the SK disk did not significantly 
tighten our constraints. Our final data vector combines 
the PIQUE Q-data and the SK T-data and thus contains 
144 + 288 = 432 pixels. 



B. Method 

W e co mpute the six power spectra describe in Sec- 
tio n jll A| using quadratic estimator method as described 
in [j7j, compu ting fiducial power spectra with the CMB- 
fast software p5[ using cosmological parameters from the 
concordance model from |2l|. We also perform a likeli- 
hood analysis as described below. 

A key challenge is separating the six power spectra 
(T, E, B, X, Y, Z). A generic quadratic band power esti- 
mator (a quadratic combination of T, Q and U pixels) 
will probe a weighted average of all six power spectra, so 
measurements of the six can in principle be afflicted by as 
many as (2) = 15 types of unwanted "leakage" whereby 



a measurement of one power spectrum picks up contribu- 
tions from another. In |Q it was argued that susceptibil- 
ity to systematic errors could be reduced by chosing the 
"priors" that determine the quadratic estimator method 
to have vanishing cross-polarizations, X = Y = Z = 0, 
and it was shown that this simplification came at the 
price of a very small (percent level) increase in error bars. 
In Appendix A, we show that this choice has an impor- 
tant added benefit: exploiting a parity symmetry, it elim- 
inates 14 out of the 15 leakeges, with only the much dis- 
cussed |?],|l0|,[l2] [ll| E — B leakage remaining. Note that 
whereas parity-conservation in the laws of physics pre- 
dicts Y = Z = 0, the proof given in Appendix A is purely 
geometric in nature, and holds even if parity-violating 
physics, foregrounds, etc., produce non-vanishing Y and 
Z. 

This result is useful in all but eliminating the leakage 
headache, which would otherwise complicate both calcu- 
lation and interpretation. It also has important impli- 
cations for other approaches, notably the currently pop- 
ular maximum-likelihood (ML) method as implemented 
by the MADCAP software M and applied to Maxima, 
BOOMERanG, DASI and other experiments: generaliz- 
ing it to polarization in the obvious way will produce the 
exactly the sort of leakage that our method eliminates. 
However, this problem can be eliminated as described 
below. 

The quadratic estimator (QE) method is closely re- 
lated to the ML method: the latter is simply the 
quadratic estimator method with B = F _1 in equa- 
tion (||), iterated so that the fiducial ("prior") power 
spectrum equals the measured one 28 1. For a detailed 
comparison between these two methods, see section IV. B. 
in |?]. The ML method has the advantage of not requir- 
ing any prior to be assumed. The QE method has the 
advantage of being more accurate for constraining cosmo- 
logical models — since it is quadratic rather than highly 
non-linear, the statistical properties the measured band 
power vector q can be computed analytically rather than 
approximated, which allows the likelihood function to be 
computed directly from q (as opposed to x), in terms of 
generalized x 2 -distributions [ p9| . 

Both methods are unbiased, but they may differ as 
regards error bars. The QE method can produce inaccu- 
rate error bars if the prior is inconsistent with the actual 
measurement. The ML method Fisher matrix can pro- 
duce inaccurate error bar estimates if the measured power 
spectra have substantial scatter due to noise or sample 
variance, in which case they are unlikely to describe the 
smoother true spectra. A good compromise is therefore 
to iterate the QE method once and choose the second 
prior to be a rather smooth model consistent with the 
original measurement. The lesson to take away from Ap- 
pendix A is that if the ML method is used, one should 
reset X = Y = Z = before each step in the iteration, 
thereby eliminating all leakage except between E and B. 
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Table 1 Polarization power spectrum^ 



8T 2 ± a 2 [^K 2 ] A£ ± a e ST [fiK] 



T 
E 
D 
X 
Y 



2183 ±905 
-2.4± 17.4 
-8.1± 16.0 
-145.5±222.2 
-7.9±206.2 



151.3± 85. 



46.7^ 



153.4± 87.6 < 3.9 (5.7) 
141.6± 73.7 < 2.8 (4.9) 
137.4± 69.4 < 8.6(17.3) 
151.3± 82.5 <14.1(20.1) 



(a) Results from combined PIQUE and SK data; 
'^Values in parentheses are 2-a upper limits. 



C. Results 

Table 1 shows the result of our band-power estimation. 
Here we use 50 multipole bands of width A£ = 20 for each 
of the six polarization types (T, E, B, X,Y, Z), thereby 
going out to £ — 1000, and average the measurements 
together into a single number for each polarization type 
to reduce noise. We eliminate sensitivity to offsets by 
projecting out the mean (monopole) from the T and Q 
maps separately. The values shown in parentheses in the 
most right column of Table 1 are 2-a upper limits (those 
are also the values that we present here as our upper 
limits). 

The detection of unpolarized power is seen to be con- 
sistent with that published for the full SK map || . The 
table shows that we detect no polarization or cross- 
polarization of any type, obtaining upper limits, just as 
the concordance model predicts. No results of Z are re- 
ported since PIQUE provides only Q data (note that Q 
and U are needed to isolate Z-polarization — we will 
include Z in in Section IV). 

The window functions reveal substantial leakage be- 
tween E and B, so the limits effectively constrain the 
average of these two spectra rather than both separately. 
For this reason, and to recast our constraints in terms 
of the TE correlation coefficient r^, we complement our 
band-power analysis with a likelihood analysis where we 
assume B = 0. Specifically, we set B = Y = Z = and 
take each of the remaining power spectra (T, E, X) to be 
constant out to £ = 1000. 

We first perform a simple 1-dimensional likelihood 
analysis for the parameter E using the PIQUE data alone 
(discarding the SK information), obtaining the likelihood 
function (bottom solid line in right panel of Figure |J) in 
good agreement with that published by the PIQUE team 
M. They find 95% of the area for E < 10£iK — our like- 



lihood curve drops by a factor e at a slightly lower 
value E w 7/iK as expected, since the likelihood curve is 
asymmetric and hence highly non-Gaussian. This non- 
Gaussianity also means that the precise confidence levels 
of the upper limits in Table 1 should be taken with a 
grain of salt. The B- limit is rather low because noise fluc- 
tuations give a negative best estimate (recall that what 
is measured is the sky power minus the expected noise 



power, which can be negative), and this preference for 
negative B pulls down the E estimate too since the strong 
leakage implies that it is really measuring a weighted av- 
erage of E and B. 

We then compute the likelihood function including 
both PIQUE and SK data in the 3-dimensional space 
spanned by (T,E,ri) and compute constraints on indi- 
vidual parameters or pairs by marginalizing as in 
Figure || shows that this produces a T-measurement 
T k 50//K, consistent with that for the full SK map || 
(left panel). Figure [| shows our constraints in the (E, rg)- 
plane after marginalizing over T. The "C" shape of the 
contours in here are not generic: we performed a series of 
Monte Carlo simulations, and some produced "backward 
C" shapes instead, others "S"-like shapes. This figure 
also shows that our constraints on the cross-polarization 
are weaker than the Schwarz inequality rt < 1, so in this 
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FIG. 8. Likelihood results using PIQUE Q-information alone 
(right panel, solid line) and using both PIQUE Q- and SK 
T-information and marginalizing (remaining two curves). From 
top to bottom, the two horizontal red lines correspond to 68% and 
95% of C.L., respectively. 
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FIG. 9. Joint constraints on E polarization and rg after 
marginalizing over T. From left to right, the contours show that 
the likelihood function has dropped to e~ , e~ 3 and e -4 6 times 
its maximum value, which would correspond to 68%, 95% and 99% 
limits if the likelihood were Gaussian. For comparison, the concor- 
dance model predicts (E, r^)=(0.82,-0.65) at £=137, the center of 
our window function for X (see Table 1). 
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sense the data has taught us nothing new. This funny 
situation is unique to measuring correlations, since any 
measurement of say T or E, however noisy, will always 
rule out some class of theoretically allowed models. How- 
ever, it is important to point out that this failure to beat 
the Schwartz inequality does not mean that an exper- 
iment i s far from detecting polarization: as we saw in 
Section II C , the actual correlation coefficient can be of 



order unity, so the step from overcoming the Schwartz 
bound to detecting an X-signal may in fact be quite 
small. Figure 10 summarizes all polarization limits pub- 
lished to date, and indicates that the detection of cross- 
polarization may indeed be just around the corner. 

Cross-polarization information could, in principle, re- 
duce error bars on E as well. As an extreme example, if 
we had measured that 77 w 1 on the angular scales probed 
by PIQUE, with tiny error bars, then we would know 
the exact spatial template of the .E-map from the T-map 
and could therefore fit for its amplitude quite accurately. 
The effect of adding cross-polarization information on the 
PIQUE S-limits is shown in Figure || (dashed line in right 
panel). In this case, we had no meaningful constraints 
on re, so it is not surprising that this approach does not 
help us: we see that adding the unpolarized temperature 
simply weakens the upper limit slightly courtesy of noise 
fluctuations. 

Since re is expected to oscillate between positive and 
negative values, using a flat (constant) m in the likelihood 
analysis runs the risk of failing to detect a signal that is 
actually present in the data, canceling out positive and 
negative detections at different angular scales. This is not 
likely to have been a problem in our case, since re is uni- 
formly negative in our sensitivity range £ = 137 ± 69 for 
the concordance model, but for future experiments with 
higher signal-to-noise, it will be important to parametrize 
re in a more physical way - either with separate bandpow- 
ers in multiple bands or directly in terms of cosmological 
parameters. 



D. Future prospects 

What improvements in data would enhance the scien- 
tific potential of an experiment like PIQUE the most? 
Reducing the noise level in the polarization measure- 
ments would obviously improve the limits on E, B and re- 
We find that substantial reductions would be needed to 
make a qualitative difference - a simulation merely dou- 
bling the amount of data, adding [/-polarization, did not 
produce an 77-detection or beat the Schwarz inequality. 
Figure |ll| shows which sort of sensitivity improvements 
are needed for a marginal detection in the E-re plane us- 
ing the same number and position of pixels as for PIQUE. 
The Schwarz bound inequality is seen to be beaten by a 
substantial margin on the high side. This illustrates the 
potential value of experimental groups coordinating to 
observe the same sky regions. Increased sky coverage 



and a more two-dimensional geometry clearly helps, and 
we explore such an example in the next section. 
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FIG. 10. Summary of upper limits on polarization so far. From 
top to bottom, the three curves show the concordance model predic- 
tions for Cj , Cf and , respectively. Four reionization models 
with t=0.1, 0.2, 0.3 and 0.4 are also plotted (left thin lines from 
bottom to top). The limits are: PW65 Q, C78 ^3), N79 Q, 
LS79 ||, LS81 §6|. 
hexagons P88 



S98 [[yii_L83 |3 i _W93 
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K01 JiH (all in upper panel) and our new upper limit on X (lower 
panel) . 




E amplitude [/-tK] 

FIG. 11. E-re likelihood for the case of lower noise (50 fiK per 
SK pixel, 5 /^K per PIQUE pixel), showing that this sensitivity 
level can give interesting constraints on r t . The data used here is 
a Monte Carlo simulation with (E, rg )=(0.82,-0.65), which is what 
the concordance model predicts at £=137, the center of our window 
function for X (see Table 1). 
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FIG. 12. Constraints in the E-rg plane from a simulation of the 
BOOMERanG experiment. From outside to inside, the likelihood 
contours are same as in Figure ^| A fiducial model with Cf = lfiK 
and C?=0 was used in this calculation. 



IV. CASE STUDY II: BOOMERanG DATA 

We begin case study II by extending the results pre- 
sented in M. In this section we quantify the ability of 
BOOMERanG to separate the Y and Z correlations us- 
ing (T, Q, U) maps. 

We pixelize our sky patch using the equal-area icosa- 
hedron method p(J at resolution levels 35, corresponding 
to 361 BOOMERanG pixels^. We apply the quadratic 
estimator described in [?j just as we did for PIQUE, with 
the fiducial power spectra Cj computed using CMBfast 
software [p5f using cosmological parameters from the con- 
cordance model from plj [. In our fiducial model we set 
Cf = l^K 2 and Cf = Cf = Cj = Cf = 0, and elim- 
inate sensitivity by offsets by projecting out the mean 
(monopole) for T, Q and U maps separately. 

Figure [l^ shows how important the unpolarized coun- 
terpart is when constraining the polarized power spec- 
trum: high signal-to-noise temperature data is seen to 
substantially improve how well E can be constrained. 
The 1/LtK E signal used in our case study is detected 
at high significance when using polarization data alone. 
Although adding T-information is obviously necessary to 
constrain re, we found that it helped only marginally for 
constraining E: we repeated our likelihood analysis us- 
ing T, Q and U data jointly, marginalizing over T and 
re, and obtained essentially unchanged error bars. 




Multipole 1 

FIG. 13. BOOMERanG power spectrum T, E, B, X, Y and Z. 
The power spectrum was calculated for 5 bands of size A£=24, Cf 
was computed using CMBfast software, and we assumed Cf=l 
and Cf = Cf = Cf = Cf = when running our Monte Carlos. 

Finally, Figure |l3| shows that we can separate the Y 
and Z correlations if we use (T, Q, U) maps jointly. We 
also verified numerically that there is no leakage between 
the TB and EB correlations, as proved in Section VI. 
Note that the results in this section underestimates the 
true power of Boomerang 2002 by only using enough pix- 
els to probe the power out to I ~ 120, the key intention 
being simply to demonstrate that the both experimental 
sensitivity and our analysis method is adequate. When 
the real data is available, a more computer intensive anal- 
ysis with or order 10 4 pixels will be worthwhile. 



V. CONCLUSIONS 

We have presented the first attempt at measuring the 
CMB cross-polarization, using the PIQUE and SK data 
sets. We obtain upper limits of TE = X < 8.6 fiK (95% 
CL) and TB = Y< 14.1 (95% CL). Limits for E are 
in concordance with the values presented in Q, and are 
much higher than the expected NCP foreground levels 
El'HO' ^ e also discuss theoretical and practical issues 



'* We use the icosahedron pixelization since it has the round- 
est (mainly hexagonal) pixels and is highly uniform. Although 
we did not use it here, the HEALPIX package is useful 
in allowing azimuthal symmetry to be exploited for saving 
computer time. 



M extrapolated an upper limit of 0.5 (iK for the polarized 
dust emission from the IRAS 100 /im map ^] and an upper 
limit of 0.4 fjK for the polarized synchrotron emission from 
the Brouw & Spoelstra Q and Haslam [ [48| maps — both 
limits were for the NCP region and the PIQUE observing 
frequency. 
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relevant to measuring cross-polarization and illustrate 
them with simulations of the upcoming BOOMERanG 
2002 experiment: we find that substantial improvements 
on measuring the polarized power spectra would be pos- 
sible if the noise could be lowered in the unpolarized map 
used for the cross-correlation. Among other things, we 
show that the well-known problem of EB leakage, which 
complicates measurements of E and B power, vanishes 
when measuring the three cross-polarization power spec- 
tra (TE, TB and EB). 

VI. APPENDIX 

In this appendix, we prove the no-leakage theorem de- 
scribed in the text. Specifically, we show that that our 
quadratic estimator method gives no leakage between any 
of the 15 power spectrum pairs except E/B. 

We follow the notation of Q throughout. We let x 
denote the vector that contains the measured tempera- 
ture and Stokes parameters at each pixel and consider 
quadratic estimators of the different power spectra, 

qi = x*Q;x = tr [Qixx*], (6) 

where i labels at the same time both the polarization 
type P (T, E, B, A, Y and Z) and the multipole I to be 
measured, qi probes a weighted average of the power 
spectra, 

(q t ) = tr [Q 4 C] 

= fo +E T,Wi, p P ,C?, (7) 

P' = l t'=2 

where b = tr [Q;£] is the contribution from experimen- 
tal noise, S is the noise covariance matrix, S is the 
covariance matrix due to the cosmological signal, and 
C = S + S is the total covariance matrix. In equation 
we introduced generalized window functions, 

Wf pi =tT[CliPvl (8) 

where P = dC/dCf. For a fixed (P,i), they show 
the expected contributions to the band power estimate 
qi not only from different £- values, but also from dif- 
ferent polarization types. For instance, an estimate of 
i?-polarization may inadvertently pick up a contribution 
from i?-polarization as well, since it is difficult to separate 
the two with only partial sky coverage ]7|, |l0| , p^ -^5| . Such 
leakage between different polarization types is clearly un- 
desirable since it complicates interpretation of the mea- 
surements. The aim of the appendix is to show that for 
our particular choice of estimator most elements of 
the window functions Wg, P , vanish, and that there is in- 
deed no other leakage than between any E and B — not, 
say, between X and Y or between E and Z. 

In j?J it was shown that the quadratic estimator defined 
by 



Qi=2M^BiiC- 1 P J -C- 1 , (9) 

j 

distills all the cosmological information from x into the 
(normally much shorter) vector q if C is the true covari- 
ance matrix. Moreover, if C is a reasonable estimate of 
the true covariance matrix, then the data compression 
step of going from x to q destroys information only to 
second order. In equation (Q), B is an arbitrary invertible 
matrix and the normalization constants A/i are chosen so 
that all window functions sum to unity. This means that 
we can interpret qi as measuring a weighted average of 
our unknown parameters, the window giving the weights. 

For our estimates we constructed the covariance matrix 
appearing in equation (^) setting Cf — Cj = Cf = 
and Cf = Cf . As explained in detail in 0, the 
choice Cf = Cf = was made because it makes 
the estimates of the temperature power spectra only 
depend on quadratic combinations involving two mea- 
sured temperatures, the estimates of X and Y only de- 
pend on quadratic combinations involving one Stokes pa- 
rameter and one temperature and the estimate of the 
E, B and Z spectra only depend on products of two 
Stokes parameters (otherwise additional non-intuitive 
terms get included, say temperature autocorrelations 
when measuring E, increasing susceptibility to system- 
atic errors). These facts immediately imply that the only 
mixed windows that could potentially be non-zero are 
Wf%, Wf, E z , Wf%, i.e., that of the ( 6 2 ) = 15 types 
of potential leakage, only EB, XY, EZ and BZ leakage is 
possible. What we will show in this appendix using an 
argument based on the parity of the different fields is 
that only W^, E B is non-zero for our choice of C. 

We consider a parity transformation such that the co- 
ordinate r' = — r, where primes indicate the coordinates 
after the transformation. We define the operator II to 
be the one that transforms the vector of temperatures 
and Stokes parameters x under the parity transforma- 
tion, x' = IIx. It satisfies II 2 = I, the identity. For 
example, if the Stokes parameters at each point on the 
sky were defined with respect to the spherical coordi- 
nate system, then II transforms them as Q' = Q and 
U' = —V. Under a parity transformation, we also have 
V = T, E' = E and B' = -B. 

To understand how the matrices C and P behave un- 
der parity, we need to understand their structure. Let us 
consider an arbitrary pair of points labeled A and B. The 
matrices are most easily described when the Stokes pa- 
rameters are measured with respect to the natural frame, 
that is when Q is defined as the differences in intensity 
in the directions parallel and perpendicular with respect 
to the great circle that joins A and B. In that case the 
structure of the covariance matrix is, 

(T A T B ) = J2 f(Cf 
i 

{Q A Q B )=J2f2Cf + ^Cf 
i 



(0) 
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(UaU b ) 
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{TaQb) 




(T a Ub) 




{UaQb) 


= E(/s-/3)Q z 



(10) 
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where the coefficients ff are known functions of the an- 
gular separation between A and B. From the above ex- 
pressions we can calculate what the different P-matrices 
are in this frame. The matrix P^ has non-zero entries 
only for the terms involving Ta and Tb , Pe and Pb only 
terms involving Qa Qb and Ua Ub, Px only for those 
involving Qa Tb and Qb Ta, Py only for those involving 
Ua Tb and Ub Ta and Pz only for those involving Qa 
Ub and Qb Ua- Moreover if when constructing C we 
took Cf and Cj to be zero, C is block diagonal, with 
no terms that mix T and Q or T and U . It is clear then 
why our estimators of X and Y only contain TQ and TU 
terms. 

Under a parity transformation in this coordinate sys- 
tem, the Stokes parameters transform as T 1 — T , Q' = Q 
and U' = —U. Under parity in (JlCf), the last two equa- 
tions change sign while the others remain the same. Thus 
if C was constructed assuming that all the cross spectra 
were zero, which implies that the last three expectation 
values in equation ([l0]) are zero, then PECII = C which 
implies that IIC _1 n = C _1 . Furthermore because only 
the last two equations in (|l0|) change sign, the P matri- 
ces satisfy np T n = p T , iip^n = p^, up^n = p B , 

nP x n = P x , nPyll = -P Y and UP Z U = P z In 

other words, the Y and Z spectra are odd under parity 
while the T, E and X ones are even. 

We now have all the necessary ingredients to show that 
any cross window function that involves spectra with dif- 
ferent parity will be zero. These window functions will 
be given by the values of tr [C _1 PjC _1 Pi/]. If the two 
P-matrices have different parities, we have 

tr [c^np.nc^np^n] = -tr [c^PiCr^]. (11) 

Using the fact that C _1 commutes with II and the cyclic 
property of the trace, we get 



tr [C-ip.C-^] = -tr [CT^e-^v]. 



(12) 



Thus if Pi and Pj/ have different parities, the trace will 
be zero. With our choice of C, the window functions 
mixing spectra with different parities are identically zero. 
This means that there is no leakage between X and Y, 
between E and Z or between B and Z, so the only non- 
zero mixed window function is Wf,g. In other words, 
out of the 15 potential leakages, our method eliminates 
all except that between E and B. 
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